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We consider a Fermion in the presence of a rotating black hole immersed in a universe with 
positive cosmological constant. After deriving new formulae for the event, Cauchy and cosmological 
horizons we adopt the Carter tetrad to separate the aforementioned equation into a radial and 
angular equation. We show how the Chandrasekhar ansatz leads to the construction of a symmetry 
operator that can be interpreted as the square root of the squared total angular momentum operator. 
Furthermore, we prove that the the spectrum of the angular operator is discrete and consists of simple 
eigenvalues and by means of the functional Bethe ansatz method we also derive a set of necessary 
and sufficient conditions for the angular operator to have polynomial solutions. Finally, we show 
that there exist no bound states for the Dirac equation in the non-extreme case. 

PACS numbers: Valid PACS appear here 


I. INTRODUCTION 

In this paper we study the spectral properties of the angular operator associated to massive Dirac particles outside 
the event horizon of a non-extreme Kerr-de Sitter (KdS) manifold and we prove the absence of bound states. The 
KdS metric is a solution of Einstein field equations describing an asymptotically de Sitter space-time containing a 
rotating black hole [Ij . Although it is not the most general model of the exterior region of a black hole we can analyze 
theoretically, since the black hole charge is not taken into account, it represents indeed the most realistic model in 
astrophysics because in general black holes are embedded in environments filled with gas and plasma and, hence any 
net charge is neutralized by the ambient matter. Moreover, the Wilkinson Microwave Anisotropy Probe (WMAP) 
indicated that our universe contains a dark energy component equivalent to a tiny positive cosmological constant [H- 
Hence, it is more than reasonable to study the Dirac equation in the geometry of a non-extreme KdS black hole. A 
first attempt to analyze the Dirac equation in the aforementioned metric traces back to where the separation of 
the Dirac equation into an angular and a radial system was achieved by using the Kinnersley tetrad coupled with 
a Chandrasekhar-like ansatz for the spinor. However, no physical explanation was given concerning the symmetry 
operator lurking behind this so fruitful ansatz. We fill this gap by showing how the Chandrasekhar ansatz leads to the 
construction of such an operator that can be interpreted as the square root of the squared total angular momentum 
operator. Separation of the Dirac equation into ordinary differential equations in general higher dimensional Kerr- 
NUT-de Sitter space-times and in type D metric has been studied in Q and 0 , respectively. Furthermore, @ proved 
the essential self-adjointness of the angular operator and the absence of normalizable time-periodic solutions for the 
Dirac equation in the KdS metric. In that regard we extend the results of @ by proving the self-adjointness of the 
angular operator. Furthermore, we give a rigorous derivation of the spectrum of the angular operator and also derive 
a set of necessary and sufficient conditions so that the angular system admits polynomial solutions. Last but not least, 
we offer a proof of the absence of bound states which is shorter and relies on a different method than that adopted 
by 0. 

The paper is organized as follows. In Section 1 we gwe a short introduction motivating the importance of our findings. 
In Section 2 we complement the results given in [ij by giving a thorough classification of the roots of the quartic 
polynomial equation controlling the location of the horizons paying special attention to their algebraic multiplicities. 
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In the non-extreme case of a KdS black hole we construct expansions for the positions of the horizons with respect 
to the small parameter A representing the cosmological constant. We also show that there are two different types 
of extreme KdS black holes and for each of them we obtain analytical formulae for the horizons. We conclude this 
section by analyzing the naked case which is characterized by fixed values of the rotation parameter and the mass of 
the black hole. Also in this scenario we are able to give the exact position of the cosmological horizon. In Section 3 
we employ Carter tetrad and a Chandrasekhar-like ansatz to separate the Dirac equation into an angular and radial 
system. The choice of the Carter tetrad is motivated by the fact it allows to give a more elegant treatment of the 
separation problem and at the same time it leads to a simpler form for the radial and angular systems than those 
derived by who instead used the Kinnersley tetrad. We also construct a symmetry operator of the Dirac equation 
in the KdS metric generalizing the one obtained in [13 and show that it can be interpreted as the square root of the 
squared total angular momentum operator. In Section 4 we study the angular eigenvalue problem. More precisely, 
we prove the self-adjointness of the angular operator and we show that its spectrum is purely discrete and simple 
by using an off-diagonalization method. Furthermore, we derive a set of necessary and sufficient conditions for the 
existence of polynomial solutions of the angular system. Finally, in Section 5 we show that the Dirac equation in the 
KdS metric does not admit bound states solutions. 


II. THE KERR-DE SITTER METRIC 


The Kerr-de Sitter (KdS) metric represents a rotating black hole immersed in an asymptotically de Sitter space-time 
with positive cosmological constant A. In Boyer-Lindquist coordinates (t, r, d, (p) with r>0, 0 <i 5 < 7 r, 0<(/3<27r 
this metric is given by SHI [13 

, 9 / . . 9 „ , n9 ^ ^dr^ d{P\ A,5Bing'd. . / 9 9n , ,9 /.n 

ds^ = — asm^ "ddip) — S ( — -h — — I -— [adt — (r^ -|- a?)dip\ (1) 

V Af. A.^ J 

with 

= {r^ + a^) — 2Mr, A,j = 1cos^ S = l-|-^^^^, S = -|-cos^ 

\ o J o o 


where M and a are the mass and the angular momentum per unit mass of the black hole, respectively. By 
d^Kdsi-^, M,a) denote the family of Kerr-de Sitter space-times. Observe that if A = 0 the above metric re¬ 
duces to the usual Kerr metric. For A 0 and a = 0 we obtain the line element of a Schwarzschild-de Sitter black 
hole. More over, if A ^ 0 and M = 0 = a the manifold becomes a de Sitter universe with cosmological horizon at 
Tc = -\/3/A. Finally, if M = 0 and a ^ 0 7 ^ A it has been shown by [l| that the coordinate transformation 


T = t/2 


_ aA 


r cosi9 = 2/COS0, = 


a? sin^ -d 


maps the family of space-times FKds{A,Q,a) to a family of de Sitter universes. Since A > 0 it follows that is 
always positive and the position of the horizons is determined by the roots of the quartic equation A,. = 0. We will 
use the complete root classification method developed by in order to study the zeros of this equation. To this 

purpose we rewrite it as 


+ pr^ -I- gr -1- u = 0 


( 2 ) 


with 



q 


6 M 



( 3 ) 


According to 11^ e quation ([2]) will have four, two, or no real roots. However, a more subtle root classification than 
that offered by [1§ | will arise because of the different algebraic multiplicities of these roots. Furthermore, we recall 
that in 1998 [T^ and [13 used Type la supernovae to show that the universe is accelerating, thus providing the first 
direct evidence that A is non-zero, with A Ri 1.7 x 10“^^^ Planck units. We will further consider A as fixed. The 
analysis of the roots of equation Q is greatly simplified if we rescale the radial variable according to p = rjrc,dS 
where rc,dS = ■\/3/A is the cosmological horizon of the corresponding de Sitter universe belonging to the family of 
space-times J^Kds{A,0,0). Then, our equation ([5]) becomes 


p^ -I- pp^ -I- gp -I- M = 0 


( 4 ) 
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FIG. 1: Allowed values for the parameters m and x = 2a^. The brighter shaded region corresponds to the inequality L > 0 
whereas the darker shaded region to <5 > 0. This demonstrates that 5 > 0 is the more stringent inequality. 


with 


p = —(1 —a^), <7 =/i, u =—a^ 


where a = ajrc^dS and ^ = 2M/rc,dS- Taking into account that the Schwarzschild radius = 2M is always smaller 
than the de Sitter cosmological horizon we find that the parameter /r can vary only on the interval [0,1). As in MM 
we introduce the auxiliary polynomials 

S{p, q, u) = 256tt^ — I28j^u^ + 144pg^u + IQp'^u — 27'(f — 4p^g^, 

L(p, q, u) = 8^ — — 2p^. 

Then, we have the following root classification for the quartic equation (|T|. 

• There will be four real distinct roots if <5 > 0, and L > 0. These two conditions give rise to the following system 
of inequalities 

27 //^ 

— —+ (1 + 33q;^ — 33q;^ — ol ^)^— — Aa^ — 6a® — 4a® — a^® > 0, (5) 

+ 2 + 2a^ - 2a^ - 2a® > 0. (6) 

The first inequality in ([5]) can be cast into the form 

0 < fm^ + g{a)m + h{a) =: P 2 {rn) (7) 

with / = —9, g{a) = (1 — a^)(a"‘ + 34a^ + 1), and h{a) = —3a^(l + a^)'*. Note that ([7]) is a polynomial 
inequality of degree two in the parameter m := M^A. The inequality |6] takes the form 

6m < 1 + a^ — a'^ — a® =: Ps^a) (8) 

A straightforward analysis of these inequalities is presented in Figured] which demonstrates the allowed region 
in the variables a and m according to © and ([5|). It is obvious that there exist extremal values for a and m 
(which we will discuss below). In the same figure we see that the allowed region of the inequality (5 > 0 lies 
inside the allowed region of the inequality L > 0. In principle, it would suffice to consider only (I7|). It is also 
possible to obtain some analytical results. First of all, we observe that P 2 {m) is a concave down parabola with 
respect to the parameter m. The solution set of the inequality P 2 (jn) > 0 will be non empty if the polynomial 
P 2 has two zeros. The condition for that reads 

g^ + A\f\h = (a^ — 4a + l)®(a^ + 4a + 1)® > 0 (9) 


which will be satisfied for a < —2 — or — 2 + < a < 2 — or a > 2 + and the zeroes are given by 

—g[a) ± •yg^ap”+4)7]7i(^ —(1 — a^)(a^ + 34a^ + 1) ± (a^ — 4a + l)^(a^ + 4a + 1)^ 


2 / 


• ( 10 ) 


m±{a) 


18 











4 



FIG. 2: Plot of the roots of the polynomial P 2 on the interval a € (—2 + %/3, 2 — v^). The solid line corresponds to m_ while 
the dotted line represents m+. 


It is not difficult to verify that m+{a) is negative if a < —2 — \/Z or a > 2 + whereas it is positive 
whenever —2 + < a < 2 — Since m+(a) is even we have that on the latter interval 0 < m+(a) < 

416/\/3 — 240 = 0.177709. Clearly, m_(a) is also even and positive on the interval —2 + < a < 2 — \/3 

where 1/9 < m-{a) < 416/\/3 — 240. From Fig. [ 5 ] we see that m+(a) < m < m-{a). The maximal value of 
m is rumax = 416/\/3 — 240 = 0.177709 and is attained at a = —2 + \/3 and a = 2 — v^. Furthermore, the 
minimal value of m is rrimin = 1/9 = 0.111111 and is realized for a = 0. To summarize, from the inequalities 
© and ([9]) we derive absolute limits on a and M, more precisely 


and 


O'min — O^max ^ CL ^ Clmax 


2\/3-3 


0.464101 

^7^ 


Mr, 


1 

sVa 


<M < 


vx 



( 11 ) 


( 12 ) 


Of course, given an a or a the mass can range between rrimax and rrimin and viceversa, for a given M or 

m we get a range for a or a. Descartes’ rule of signs predicts that m will have one negative root at r _ 

without physical meaning and three positive roots at r_, r+, and Tc representing the Cauchy, the event, and 

the cosmological horizon, respectively. Note that r _< 0 < r_ < r+ < rc. Moreover, we have the factorization 

Ar = —(r — r _ ){r — r_)(r — 'r+)(r — r,.) and > 0 on the intervals 0 < r < r_ and r_|_ < r < Tc- Note 

that the presence of r _can always be removed since we can express this root as r_= — (r_ + r_|_ + r^) by 

using one of Vieta’s formulae. For completeness we give also analytical expressions for the roots of m which 
are represented by 




- 2:1 


y/-Z2 + \/-^3 


--21 


'-^1 - 'J-Z-i 


r+ = 


'-z\ - J-Z 2 


^-z-i 


Ta = 


~\/—Z\ + \/ — Z2 + 'J—Z'i 


(13) 


where Zi = Wi + (2p)/3 for every i = 1, 2, 3 and 


Wl = + ^2, 11)2 = UJ + UJ^2, W3 = Ul^i + CO ^2 


(14) 
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with w = exp(27ri/3) and 



The details of the derivation of the above formulae can be found in the appendix. Note that the formulae for the 
roots are too unwieldy for general use. However, since A can be interpreted as a small parameter, we can use 
perturbation theory of algebraic equations to derive simpler formulae for the Cauchy, event, and cosmological 
horizons. To this purpose, let us consider the polynomial 

Ka? \ n n 

1-^ j + 2Mr — a?. 

Then, we immediately see that it can be written as an unperturbed polynomial plus a perturbation containing 
A, i.e. Pa(?’) = 7b(?') + E\{r) with To(r) = —+ 2Mr — a^, and E\{r) = Ar^(r^ + a^)/3. Clearly, the roots 
of To(r) = 0 are given by r±^K = M ± Vwhere r+^K and r-^K coincide with the event and Cauchy 
horizon of a Kerr black hole, respectively. Since E\{r) vanishes as A tends to zero, a lemma in perturbation 
theory (see page 83 in [13) ensures that the equation Pa(?') = 0 has at least two roots r_(A) and r+(A) such 
that r±(A) —7> r±^K as A —>• 0. Substituting r±_A = ?'±,a: + ?’i,±A + 0{A?) into the equation Pa(?’) = 0 yields 

To{r±^K) + ^ + a^) - 6(r±,if - Af)ri,±] A + ^(A^) = 0. 

From the definition of r±^K it follows immediately that Tq(j'±^k) = 0. Moreover, the fundamental theorem of 
perturbation theory implies that the coefficients of the various powers of A can be set to zero. Hence, we find 
that the coefficients ri^± are given by 


PA{r) = -Ar = - 


A 

"s' 


“ 6{r±^K -M) ~ 3VM2 - a? ' 

Note that for A approaching zero, inequality ([7]) reduces to the condition of a non extreme Kerr black hole, that 
is > a^, and therefore the coefficients ri^± are always well defined. Hence, the event and Cauchy horizons 
of a non extreme Kerr-de Sitter black hole can be represented as follows 


r± = r±^K ± 


Mrlj, 

3VM2 - a2 


A + C>(A2). 


Concerning the cosmological horizon we rescale the radial variable as r = w/\/~A and we consider the equation 

Pa(A“^/^w) = — A~^vP' + 2MK~^/‘^w — = 0. 

o o 

Multiplying throughout by A and setting /3 = A^/^ we obtain the regular perturbation problem 

- Zv? + QMPw - = 0. (15) 

As /3 —>■ 0 we get the unperturbed equation = 0 having roots at 0 and ±-\/3. Taking into account that 

for vanishing a and M the Kerr-de Sitter metric represents a de Sitter universe with cosmological horizon at 
rc,dS = \/3/A, we will look at an expansion of the form w = -I- wi/3 -I- ^ 2 /?^ + Substituting 

such an expansion into (fTKI) and setting the coefficients of the various powers of f3 to zero, we obtain wi = —M, 
W 2 = —V^M'^j2, and W 3 = —M(M^ — a^)/3. Finally, we conclude that the cosmological horizon for a non 
extreme Kerr-de Sitter black hole admits the following expansion in A, namely 

re = re.ds - mVA - ^^A - y (M^ - + 0{A% 


In the case a —>■ 0 the above expression gives the position of the cosmological horizon for a Schwarzschild-de 
Sitter black hole. 
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• There will be two distinct and two coinciding roots whenever L > 0, <5 = 0 and —amax < a < amax- The 
value of m can be only on one of the two branches m± (a). Since the condition (5 = 0 reduces for A —>■ 0 to the 
constraint = a^, we will call this case the extreme Kerr-de Sitter black hole. Let a G (—2 + \/3, 2 — \/3). 
Since ([S]) is biquadratic in the corresponding roots are 

Aii,±(a) = F{a) + G{a), /i2,±(a) = - ^\/G{a) 

with F{a) = &g{a), G{a) = (a^ — and a± = 2 ± '/S. We have the following two cases 

— if ^ the event and cosmological horizons coincide, i.e. p+ = Pc, and 


= i V 6 — 6 q;^ + 6 \/— 14a^ + 1, 


P- 


This implies that 


-pc + \/l - - 2 pl. 



^ ^6 - 2a2A + 2Va4A2 -42a2A + 9 


r_ 



with expansions in the small parameter A given by 

r, = ^ r_ = + oiA^/^). 

A 3 9 2 8 


If ^ = p 2 ,+ io:), the Cauchy and event horizons coincide, i.e. p- = p+, and 


p+ = - Y ^6 — 60(2 _ 6 -\/a^ — 140(2 +T, pc = —p+ + \ j I — — 2 p \, 


that is 


1 /3 


r+ = 


6 - 2a?A - 2\/a^A^ -A2a?A + 9 


6 V A L 

with expansions in the small parameter A given by 




rc = -r+ + \l - a? - 2r\ 


41^ 


r+=a + -a^A + C>(A2), = \ — - a - —a^\fA- -a^A -_p 0{A^). 


72 


For the details of the derivation we refer to the appendix. From the results above we see that the extreme 
Kerr-de Sitter black hole in reality contains two possible scenarios, li p = pi^+ we have < 0 in the region 
between the Cauchy and the cosmological horizons, thus signalizing the presence of a dynamical universe where 
the central singularity is shielded by a Cauchy horizon. A similar scenario emerged also from the analysis of 
0 where the horizons of Reissner-Nordstrom-de Sitter black holes have been investigated. In the limit a —>■ 0 
we have M = 1/(3-\/A) and this case corresponds to an extreme Schwarzschild-de Sitter black hole where the 
event and cosmological horizons coincide. In this case we have a dynamical universe exhibiting a naked central 
singularity. Moreover, if ^ = p 2 ,+, we have instead > 0 in the region between the event and cosmological 
horizons. This case reduces to the usual extreme Kerr case in the limit A —?> 0 as we can see from the above 
expansion for r^_. Last but not least, it is pretty amazing that we can find relatively simple formulae for the 
positions of the horizons in the case of an extreme Kerr-de Sitter black hole. 

• There will be two real roots with algebraic multiplicities one and three, respectively, whenever [ij, Hi 

5 = 0, L = 0, p<0, qy^O. (16) 

For L = 0 and 5 = 0 we have a = otmax and m = rrimax- The condition q ^ Q translates simply into the 
requirement M ^ 0 whereas the condition p < 0 is satisfied for Ao^/3 < 1. Moreover, we see that 5 = 0 and 
L = 0 at the points where the graph of p+ intersects the graphs of and P 2 ,+ - This happens for a = ±Q!_ 
with Of- = 2 — for which 


p±(±a-) = ^i,+ (±a_) = /i2,+ (±a-) = ^\/26v^- 45, 
























or equivalently, in terms of the parameters a and M 
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2\/3-3 4 / 78 V 3 -135 

Furthermore, the roots of (HI) are p _= — \/6\73^^ and p- = p+ = pc = (2/3)-\/3 — 1 to which correspond 

the following roots of namely 


18\/3-27 /2\/3- 3 

—-— , = = 

Taking into account that r > 0, it is straightforward to verify that A^. > 0 for 0 < r < Tc. Hence, this 
scenario corresponds to a naked singularity at r = 0 in the equatorial plane [Tsj l immersed in an universe with 
cosmological horizon. It is interesting to observe that this case does not admit a Kerr-like counterpart in the 
limit A —^ 0 because in the limit of a vanishing constant the condition L — Q gives rise to the contradiction 
1 = 0 . 

• The case of two real roots both having algebraic multiplicity two which corresponds to the set of conditions 
S = L — q = 0, and p < 0 does not represent a black hole since the condition q = 0 requires that M = 0 and 
therefore it will not be further analyzed. The same conclusion holds for the case d < 0 corresponding to the 
scenario of two complex conjugate roots and two real roots each having algebraic multiplicity one. 

• liS = L= p = 0 there is only one root with algebraic multiplicity four. The condition p = 0 implies that (P|) 

becomes + qr + u = 0 with q = 6M/A and u = —9/A^. On the other side ([2]) could be written as (r —ro)^ = 0. 
A quick comparison with r'^ + qr + u = 0 shows that u = = 0 which is impossible since u = —9/A^ ^ 0. 

Hence, this case can never occur. 

• If (5 = 0 and L < 0 there is only one root with algebraic multiplicity two and two complex conjugate roots. This 

means that ([2]) must be of the form (r — ro)^(r^ -I- = 0 with > 0. Comparing the latter with ([^l yields 

ro = 0, A^ = p, q = 0, and u = 0. The last two conditions imply M = a = 0. Hence, p = —3/A and since we 
are interested in the case of a positive cosmological constant, this would require that < 0. Therefore, this 
case might be relevant in the case of an anti-de Sitter space and it will not be further investigated here. 

• There will be no real roots when [d > 0 and (L < 0 or p > 0)] or [d = 0, L = 0 and p > 0]. The absence of 
real roots implies that there will be no cosmological horizon and therefore no de Sitter space asymptotically 
at infinity. However, this will never happen. First, d > 0 and p > 0 are not compatible. From the first we 
obtain c? < 0.0718 (see the analysis above) whereas the second inequality is equivalent to c? > 1. It remains 
to look into the <5 > 0 and L <0 case. Numerically the second inequality is satisfied for > 41/3 which again 
contradicts d > 0. 

Finally, after the analysis of the horizons some general comments are in order. Let us note that it is not always 
obvious that choosing one parameter in the theory constraints the choice of other parameters. In the Kerr-de Sitter 
black hole this interplay happens for a and M. Moreover, equations (ED and m sets absolute limits on a and M. 
This seems to be a general characteristic of theories containing the cosmological constant A. Indeed, several cases can 
be quoted 

1. The general condition for the existence of horizons in the Schwarzschild-de Sitter metric is [ 2 ^ 

> M (17) 

3VA “ 

2. The condition for the Newtonian limit to exist in the Schwarzschild-de Sitter case is [2l| 

» M (18) 




3. The effective potential which determines the orbits in the Schwarzschild-de Sitter case has three local extrema. 
To avoid that the first maximum coincides with the minimum one has to respect that the angular momentum 
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per mass I is smaller than •\/3?’s with Vg being the Schwarzschild radius. On the other hand insisting that the 
minimum and the second maximum do not degenerate we have 


I < 




(19) 


There exists a maximum radius of the order (the location of the second maximum) beyond which 

bound states are not possible. 

4. In the Schwarzschild-de Sitter metric different concepts of equilibrium demand a constraint on mass or density. 
The hydrostatic equilibrium demands that 

> M (20) 

9VA ~ 

The Tolman-Oppenheimer-Volkoff equation states that p^ 

A < 47rp, (21) 

with p the average density. On the other hand the virial equilibrium applies if p3 | 

A 

P ^ -^Pvacj Pvac — "Z {^^) 

on 

with p being the density and A an expression which depends on the geometry of the object. 


5. Probing the black hole evaporation via a generalized uncertainty principle one finds that there exists a maximum 
temperature corresponding to a minimum black hole mass (black hole remnant) pflj | and, in case that the 
cosmological constant is non-zero, that we have a minimum temperature corresponding to a maximum mass, 
i.e., 


T 

max 


rripi 

T 

min 

"n/A -O- Mrnax 

m-li 

^/A 


where we have restored the Planck mass nipi. 


(23) 


III. THE DIRAC EQUATION IN THE KERR-DE SITTER METRIC 


A fermion of mass rUe and charge e is described by the Dirac equation (h = c = G = 1) 


V^cI>a = -^Xa'. ^AXA' = -k<t>A 






(24) 


where denotes covariant differentiation, and {4>^,x^ ) s-re the two-component spinors representing the wave 

function. According to p^, at each point of the space-time we can associate to the spinor basis a null tetrad 
(1, n, m, m) obeying the normalization and orthogonality relations 


l n = l, m m = —1, 1 • m = 1 • m = n • m = n • m = 0. 


(25) 


Moreover, to any tetrad we can associate a unitary spin-frame (o^, defined uniquely up to an overall sign factor 
by the relations = n“, = m“, = m“, and = 1 p^. As in [l^ we denote by 0o 

and 4 >i the components of in the spin-frame (o"^, i^), and by Xo' Xi' the components of x^ in (b'^U'^)) more 

precisely po = (pAO"^, 4>i = Xo' = Xa'^^ ) ^nd Xi' = Xa''^^ i ^-nd we introduce the spinor pj = (Fi, F 2 , Gi, G 2 )"’" 

with components Fi = —po, F 2 = (pi, Gi = iXi', and G 2 = iXo'■ Then, the two equations in (1241) can be written in 
matrix form as 


/ -rUe 

0 

a+ 


\ 

0 

-rUe 


a- 

3_ 

-/3+ 

-TTle 

0 



a+ 

0 

-ITle 

/ 


= 0 , Od 


(26) 
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with 


a+ = —iV2{D + e — p), /3^ = iV2(6 + tt — a), f3-= i\/2{S + P — t), q ;_ = — i\/ 2 (A + ^ — 7 ), 

5_ = —i\/ 2 (A + ^ — 7 ), /3+= *\/2((5 +/3 — r), P-= iV2{6 + tt — a), a+=—iV2{D + e — p), 

where p, 7 , /?, r, p, tt, a are the spin coefficients and D, A, S are the directional derivatives along the tetrad 
(1, n, m, m). In what follows we consider the Dirac equation in the presence of a non-extreme Kerr-de Sitter black 
hole. The Dirac equation in this geometry was computed and separated in Q with the help of the Kinnersley tetrad 
( 2 ^. In view of the separation of the Dirac equation we choose to work with a Carter tetrad [ 2 ^ which allows for a 
more elegant treatment of the separation problem and leads to simpler forms of the radial and angular equations than 
those derived in This symmetric null tetrad also generalizes the one employed in where the Dirac equation 
in the Kerr metric has been investigated. The line element of the Kerr-de Sitter metric given by © suggests that we 
introduce the differential forms 


K = 



dv = 

A 



^ 03 sini? 

d^, = 

/\,r. ^ 



+ a^)dip - adt], ^ 



[dt — a sin^ d dtp]. 


With the help of (5.119) in [2^ we can construct a symmetric null tetrad as follows 


Hence, we have 
' 1 


Ifj. — 


n,, = — 


ruu. = 



1 /Ar 


2 A^ 


, _K + di 

a sin^ d 

sin^ d I Ar 



0 -- 
2 A/ ’ 



ia sin d 



— 0 

2E ’ ’ V 2A^ 


S ■ ■ 

—I smu 


2 S 
2 I 


= 


d‘^ - d'^ 

—JT' 

-I- dP') 


dl + idl 


V 2 ea; ’ 

/ S(r^ -I- a^) 

i y2EA; 


V 2 


^ 0 

2 E’ ’ ysss; 





2 s’'^’ ysss; 

fSasin'i? „ j A^ 


V 2 EA 






2S ’ sin'di/2EA^ J 




It is not difficult to verify that this tetrad satisfies the conditions in (|25|) and it is made of null vectors, i.e. = 
= rriqiin^ = rfifjrn^ = 0. Using the above tetrad and (2.3a) — (2.3d) in 0 the spin coefficients are found to be 
K = a = X = v = 0, p = p, ^ = 6, P = —a, TT = —T and 


d = 



^ 2 i ^ 2 V 2 EA; 


a = 


A^ / A^ ( „ fa sin 

— . + \\ —TT cot w -I-^— 

\ 2T, \ p 


fasin'i? A^ 


P 


2S’ 


where p = r -|- fa cos'd, prime denotes differentiation with respect to the radial variable and dot means differentiation 
with respect to the angular variable d. We are now ready to give a more explicit form to the Dirac equation (I26I) . If 
we introduce the following operators 


^ d E. 


the entries of the matrix in 


d 


d 


are computed to be 


. d 1 i:i ( . d d 

U± = + 7: cot d =F ^ a sm d— -f esc'd— 

ad 2 A^ \ at ap 


a± = + f{r,d)), p± =+ g{r,d)), 

FAr 


a± = ±i\l — {V± + f{r,d)), P± =-i\J—{C^ + g{r,d)), 




g(r,d) = 


1 j A^ fa sin d 

2 1 p 


with 
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As in M , we replace (1261) by a modified but equivalent equation 

Wipit,r,'d,ip) =0, W = TS-^OdS, (27) 

where ■0 = S'”= {Fi, F 2 ,Gi,G 2 )'^ and F and S are non singular 4x4 matrices, whose elements may depend on 
the variables r and Proceeding as in Lemma 2.1 in m it can be showed that for r+ < r < Tq there exist non 
singular 4x4 matrices 

S = (A^Ai,)"^/‘‘diag(p“^/^,p“^/^,p F = diag(-ip, ip, ip,-jp) (28) 

with det(S) = (SA^A,?)”^ and det(F) = such that the operator W decomposes into the sum of an operator 
containing only derivatives respect to the variables t, r and tp and of an operator involving only derivatives respect 

w = (29) 


iWer 0 yfK^Vj^ 0 \ 

0 —imeT 0 v^ArP- 

0 —imeT 0 ’ 

0 y^ArV^ 0 irUer J 

aTOeCosd 0 0 —\/A^C- \ 

0 —anie cos'd ^/A^£+ 0 

0 arrie cos'd 0 

^/A^C^ 0 0 —arrieCOS'd J 


to t, 'd and (p. More precisely, we have 

with 

/ 

V 
/ 

V 


Remark III.l It should be 'noted that this decomposition continues to hold even for extreme Kerr-de Sitter black 
holes. In addition, since the Kerr - de Sitter metric goes over into the de Sitter metric for M = a = Q, we obtain 
from the decomposition above a similar result for the Dirac equation in a de Sitter universe. Taking into account that 
in this case A,j = 1 and Ar = r^(l — Ar"^/3) with r G ( 0 , ^J‘i/A) we can find non singular 4x4 matrices 


S = diag{r ^/^), T = diag{—ir,ir,ir,—ir) 


with det{S) = (r^A^.) ^ and detfT) = such that the operator decomposes as follows 






w, 


idS) 


imer 0 Y ArT>^ 0 

0 —imer 0 \J A^D- 

0 —imeV 0 

y 0 Y^ ArT>+ 0 imeV 


' ArV_ 


W, 


(dS) 


/ 0 0 0 \ 

0 0 £+ 0 

0 0 0 

V £ + 0 0 0 J 


where 


d 1 d ^ d 1 „ B 


We compute now the commutation relations needed to construct a symmetry operator of the Dirac equation in the 
Kerr-de Sitter metric generalizing the one obtained in [l^ for the same equation in the Kerr metric. To this purpose, 
let ^ (fl) with n = R X (r+, Tc) x [0, tt) x [0, 27 r]. Then, it is not difficult to verify that the following commutators 

hold 


[W(t.r-.^), = 0 , [W(t,,,^), W] = W] = 0 . 
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Moreover, the matrix F splits into the sum F = F(r) + F(, 3 ) with F(r) = idiag(—r, r, r, —r) and F(, 9 ) = 
a diag(— cosd, cos'd, — cos'd, cos'd) satisfying the commutation relations 

[r(r-),r(,3)] =0, = [F(,j),>V(t_r.v3)] =0. 

Since the Kerr-de Sitter metric is axially symmetric, it is natural to make the following ansatz for the spinors ip 
entering in (j27p . namely 


^{t, r, <p) = 


(30) 


where ui and fc G Z are the energy and the azimuthal quantum number of the particle,respectively, and ip{r,'d) G C^. 
Inserting (l30l) in (l29l) . it can be verified that ip{r,d) satisfies the equation 


where 


with 


where 


W(,) = 


>V(^) = 


(W(,) + W(,,)) ip = o. 


( imeT 0 0 \ 

0 —iuier 0 V ArV- 

0 —iuier 0 

V 0 ^ ArV^ 0 irrier / 

( arae cos d 0 0 ^ — V AsC- \ 

0 —ame cos'd \fA^C+ 0 

0 —V A^C- arrie cos 'd 0 

\ \fAaC+ 0 0 — ame cos'd / 


dr ^^Ar . 


d 

uj(r^ + a^) + ak , C± = + A±(d) 


(31) 


(32) 


(33) 


A± (d) = - cot d ± WjJ sin 'd - 
2 A^ \ 


sind 


fc = fc + 


At this point it is instructive to compare the above expressions with the operators F>o, 2^i/2i and represented 
by equations (2.6) — (2.9) in 32| where the separation of the Dirac equation in the Kerr-de Sitter metric was achieved 
by using the Kinnersley tetrad. It can be immediately observed that the major benefit of using the Carter tetrad is 
to transform away the terms 


1 dAr 


2 Aj. dr ’ 2 ^/Ag sin 'd d'd 


(\/A)j'sind) 


appearing in (2.6) — (2.9) in [^. This in turn will lead to a very simplified form of the radial and angular systems. 
Let us set 


iP{r,d) = {R.{r)S+{d),R+{r)S.{d),R+{r)S+{d),R.{r)S.{d))'^. (34) 

According to Chandrasekhar Ansatz and proceeding as in 0 equation m splits into the following two systems 
of linear first order differential equations, namely 


/ y/APV_ -im^r- \ ^ ^ _ q 

\ imer — A ^/AP-Vj,. ) \ R+ j 


(35) 


/ —y/ApC- X + ame cos'd A / “S'- A _ q 
\ A — ame cos d ^/A^C+ / \ ‘^+ / 


(36) 
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where A is a separation constant. Note that when a = A = 0 the angular eigenfunctions S± reduce to the well-known 
spin-weighted spherical harmonics whereas for a ^ 0 = A the same eigenfunctions satisfy a Heun equation [sij. 
Finally, in the general case of non vanishing values of the parameters A and a it has been shown in that S± obey 
a generalized Heun equation (GHE). Since the GHE has been scarcely studied in the mathematical literature made 
exceptions of 331 and [s^, the next section will be devoted to study the spectrum of the angular eigenvalue problem, 
the dependence of the eigenvalues upon the relevant physical parameters, and to obtain series representations for 
the eigenfunctions. We conclude this section by giving a physical interpretation to the separation constant. To this 
purpose we will use the Chandrasekhar ansatz to generate a new operator J and we will show that it commutes with 
the Dirac operator Od in the Kerr-de Sitter metric, thus being a symmetry operator for Od- It will turn out that A 
can be seen as an eigenvalue of the operator J whose interpretation will emerge from taking the limit A —>■ 0 in the 
expression for J. Since this limit coincides with the square root of the squared total angular momentum for a Dirac 
particle in the Kerr metric obtained in 0 , we can conclude that J is the squared total angular momentum for a 
Dirac particle in the Kerr-de Sitter metric. Proceeding as in M we can construct a matrix operator 


J = 


/ 0 
0 

iap COS 

V -r^^C+ 


rp^C 

iap cos 2?4 


—iapcos'd^ 

E 
0 






0 


rp^C. 

—iap cos 

0 

0 




such that J^p = \^p. Let J = SJS ^ with S defined as in (1281) . Then, 

J = 

is a symmetry operator for the formal Dirac operator Od since [Od,J] = 0 . We s kip the proof because the method 
and the computation is essentially the same as those appearing in Lemma 2.4 in [I^. Finally, letting A —?> 0 in the 
above expression reproduces the square root of the squared total angular momentum for a Dirac particle in the Kerr 
metric obtained in llOll. 


IV. THE ANGULAR EIGENVALUE PROBLEM 

We start by observing that the angular eigenvalue problem (|551) admits the discrete symmetry P : d —tt — d so 
that 

PS'+(d) = S+{tt -i9)= S-(i9), FS-(i9) = S-(7 r-i9)= S+(i9). 

This means that if we decide to eliminate 5'+ in favor of S- in (I36p to get a second order differential equation for S- 
the corresponding equation for S+ can be obtained by applying the transformation P to the equation satisfied by S -. 
It is not difficult to verify that the equation satisfied by S- is 

,/A^C+iy^C-S-) + ^ (y2 _ ^2^2 ^Qg2 ^ (37) 

A “h (xttIq cos u 


Remark IV.1 Note that in the last term of equation (3.39) in 0/ and equation (3.2) in there is a typo and the 
term A^ -I- a^rUg cos^ d should be replaced by — afrnf. cos^ d. 

The differential equation (1571) has two singularities at z? = 0, tt. An additional singularity emerges if A belongs to the 
interval {—ame, ame). In the case that A = ±ame the third singularity coincides with the singularities at tt or 0, 
respectively. In what follows we will assume that such a singularity does not belong to the interval ( 0 , 7 r). Let 



S{d) ^ 1 / 5+(z?) \ 

v/sin'd \/sind V / 


and introduce the coordinate transformation dx/did = l/i/A)?. A simple integration gives 


1 c 


dd 




0 < fc" = 


O' 1 ^2 

T-— <1, cr = -a , 
1 -I- cr 3 


(38) 
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where the integration constant has been chosen to be zero since a; = in the limit cr —>■ 0. With the help of 110.02, 
110.06, and 113.01 in we find that 




xi{'d) = -^F ^sini?, J if 0 < < 7r/2, 

~ F{cos'd,i^/a) if 7r/2 <t} <tt. 


(39) 


Here, K and F denote the complete and incomplete elliptic integral of the first kind, respectively. Furthermore, using 
110.06 in it is not difficult to verify that x(0) = 0, 


xi(7r/2) = X2(7r/2) = xo, X2(7r) = xq + xq = ' 

This means that the interval (0, tt) will be mapped by the coordinate transformation x = x(z?) to the positive interval 
/ = (0,X2(7r)). Observe that for cr —7> 0 we can use 111.02 in to show that the interval I reduces as expected to 
the interval (0,7r). Last but not least, we have 

lim Xi(i?) = d = lim X 2 {^) 


signalizing that (15^ reduces correctly to x = -d in the limit cr —>■ 0. The coordinate transformation x = x{d) can be 
inverted and expressed in terms of the Jakobi elliptic functions as follows 


i?(x) = 


Mx)=sm 1 (sn^VSx, 


= COS ^ I sn 




if 0 < X < xo, 

if Xo < X < X2. 


(40) 


Finally, the angular system can be rewritten as a Dirac system 


(W5)(x) = 


/ 0 1 '' 

.dj 

( 

1-1 0, 

dx ^ 

V V Aa; 


vWt 


a^sin^ix) + 

arrie cos-d(x) 


awsini?(x) + 

with X S / and 9{x) defined as in (|40|) . Let us rewrite the formal differential operator U as 


S{x) = XS{x) (41) 


U = 


dx vWf 


—arrie cosz9(x) 

auj sini9(x) + — * 

'' ' SI 


'd{x) 


^ - 7^ ^awsint?(x) + 
arue cosi9(x) 


(42) 


which acts on the Hilbert space L^{I,dx)'^. To simplify the following analysis we express U as the sum of the 
unbounded operator Uu and the bounded operator Ub given by 


Uu = 


0 


dx VAa; sin 'i 9 (x) 


dx V sin 19(21) 
0 


, Ub = 


0/ \ sini^fic) 

—arrie COS v(x) ’ 

arue cos t?(x) 




In the following we always assume that A, a, TOg? ^ are real and fc S Z. The minimal operator associated to li is 
_ ^^(j) with domain of definition I1(A™“) = . Let us introduce the inner product 


($, $) = / dx <i)(x)4>*(x) 


where * denotes complex conjugation followed by transposition. Then, it is not difficult to verify that A™®” is formally 
self-adjoint. Hence, by Theorem 5.4 in [s^ the operator A™®" will be closable. Let A^ denote the closure of A™®". 
We can prove the following 

Theorem IV.2 The operator A„ is self-adjoint if and only if k G ]R\(—1, 0) and a fortiori for every fc G Z. 
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Proof. We show that Uu is in the limit point case at = 0 and d = tt. A fundamental system of the differential 
equation = Q \s 

$l(x) = (^tan^^ ^ '^gV^(fc+i)ta„-i(V5^costfW) 1 ^ ^ 

— ^/a(k+^'\ta,n~^{^/acos'i)(x)) ( 0 

U 

Let us analyze the square integrability of these solutions. To this purpose let fc > 0, a G {xo,X 2 ), and /3 G (0,a;o). 
Then, we have 


4 ) 2 ( 0 ;) = (tan 


i9(a;) 


dx | 4 ) 2 (a;)p > / dx [ tan 


d(a;) 


2 /c+l 


> 


, d{x) 
dx tan —-— 


dd 


f 'i9(a) 


tan {•d/ 2 ) 
'J Atf 


= In 


2^ — 4cr cos^ {'d/2) + 2^\/:i — 2(Tsin^ d 


cos^ {d/ 2 ) 


= + 00 , 


i9(q;) 


where in the first majorization we used the fact that cosd(a;) < 0 for a; G (q!,X 2 ) together with tan(d(a;)/2) > 0 since 
0/2 G {d{a), 7 r/ 2 ) and therefore tan“^ (•ya-cosd(a;)) is negative. The second majorization is obtained by observing 
that d{a)/2 > 7r/4 implies that tan (d(a;)/2) > 0. Furthermore, d{x)/2 > 7r/4 and hence tan(d(x)/2) > 1 on the 
interval L = (d(a)/2,7r/2). This implies that (tan (■d(a;)/2))^^+^ > tan(d(a;)/2) on L. Moreover, we also have 


dx |4)i(a;)|^ > / dx I cot 


d{x) 


2 fe+l 


> 


d{x) 
dx cot —— 


L 


0 


1 d 

—= / dd cot (tl/2) = Insin — 
vS Jo 2 


m) 


= + 00 , 


where we used the fact that for d G (0, j3) we have ^J'K^ < v^. On the other hand, we also have the following 
estimates 


' dx |4>i(o;)| 2 < e^^Ofc+i) dx (^cot (^cot ^ 


2/c + l pX2 


dx < +00 


and 


[ dx | 4 > 2 (a;)P < e^^Ofc+i) f 

Jo Jo 


dx ( tan 




2 fc+l 


< ef'/^( 2 fe+i) / 


d{a) 


2 fc+l pcx 


dx < + 00 . 


This shows that in the case fc > 0 the solution $1 lies right but not left in L‘^{I,dx)‘^, whereas the solution $2 lies 
left in L^{I,dx)'^ but it does not lie right in L'^{I,dx)^. For A: < —1 the same holds true for 4)i and ^2 exchanged. 
According to Weyl’s alternative it follows that for k G M\(— 1,0) the formal differential operator Uu is in the limit 
point case both at 0 and at tt. Hence, Theorem 2.7 in [s^ ensures that the closure of A™*" is self-adjoint. To show 
that Au is not self-adjoint for k G (—1, 0) we check that $1 and $2 are in L‘^{I, dx)^, thus Uu is in the limit circle case 
both at 0 and tt. Then, again by Theorem 2.7 in [s^ the assertion follows. We give a proof for $2 G L^{I, dx)^ in the 
case fcG(—1,1/2], the remaining cases can be treated similarly. From the initial assumption we have 2fc + l G (—1,0]. 
Hence, it follows from the inequality sin ^ ^ with 7 G ( 0 , 7 r) and the monotonicity of the cosine and tangent 

functions 


dx I 4 ) 2 (x)l^ = / dx\^ 2 {x)\ 

Jo 


2 _ ^^y/a\2k+l\ 


dx tan 


d{x) 


2 fe-|-l 


= gf V^|2fc-Hl| 


dd 

^/A^ 


d^k+l 

tan — < 

2 


^^^/a\ 2 k+l\ 


dd I tan 


0 


2 fe-|-l 


= gf" V^|2fc-l-l| 


fT!-/2 


|2fc-|-l| 


cot 


2 / 


l■7r/2 


dd I cot 


d\ 


|2fe-|-l|' 


< 
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gT V^|2fc+l| 



|2fc+l| 

cot-j 





■ ^\|2fe+l| 

cot-j 


< gT>/^| 2 fc+l| 




-|2fc+l| 



< 


gT V^|2fc+l| 


^i2fe+ii r^\^^-\2k+i\ 

Jo 2 


2 



2|2fc+l| \ 

1 - |2fc+l| j 


< + 00 . 


This completes the proof. □ 


In order to find an explicit representation for the domain of the operator Au we introduce the so-called maximal 
operator associated to Uu by D(AZ^^) = {$ £ L‘^{I,dx)^ \ <& is absolutely continuous, Uu^ S L?{I,dx), := 

W„$}. Then, by Theorem 3.9 in it follows that the adjoint A* = Since for k £ IR\(—1, 0) the operator Au 

is self-adjoint by Theorem IIV.2[ we also have that Au = 

Theorem IV.3 The angular operator AS :=US withU given by and domain of definition 

D{A) = {^ £ Lf{I,dxY I S is absolutely continuous, US £ L^(/,dx)^} 

is self-adjoint if and only if k £ M\(—1,0). In this case, A is the closure of the minimal operator A"‘^^ defined by 
^ and A^^^S := US. 

Proof. Let us start by observing that D{A) = D{Aff°‘^). Let At, be the maximal operator associated with the formal 
multiplication operator Ub, that is D{Ab) = L‘^{I, dx)"^ and Ab'i’ = Ub'i’. The operator Ab is symmetric and bounded 
in the Hilbert space L‘^{I,dx)^. Hence, Theorem 4.10 (see Ch. V in [^) shows that A = Au + Ab with domain 
D{A) = D{Au) is self-adjoint if and only if Au is self-adjoint. The result follows from Theorem IIV.21 □ 

Since A is self-adjoint its spectrum a (A) must be real. 

Remark IV.4 Consider now the angular operator in the special case = 0 and introduce the formal differential 
expressions 


S = 



0 

aw sind(a;) 



V Ax 


oa;sind(x) 


0 



Then, Theorem \IV.!^ implies that for k £ ]R\(—1,0) the operator BS = fBiS with domain of definition D{B) = 
{S' £ Lfi{I,dxY I S absolutely continuous, iBS £ L^{I,dxfi} is self-adjoint and it is the closure of the minimal 
operator H™™ given by with B^™S = IBS. This implies that the operators BS 2 = ®-i-S 2 with 

D{B) = {S 2 £ L^(I,dx) I S 2 absolutely continuous, iB+S 2 £ L'^{I,dx)} and B-Si = iB_Si with D{B-) = (Si £ 

Lfi{I,dx) I Si absolutely continuous, iB_Si £ L‘^[I,dx)} are adjoint to each other so that B = 

Let us write the angular operator A as A = with D = ame cos'd{x) a bounded multiplication operator in 

L'^{I,dxfi and 


0 B 
B* 0 




aw sind(a;) 


k 

sin'd(a;) 


By using an off-dia gon alization method as in [s^ we show that A has compact resolvent. This together with Theo¬ 
rem 6.29, Ch.HI in [38| will imply that the spectrum of A consists only of isolated eigenvalues with no accumulation 
points (— 00 , 00 ). First, we prove that the discrete spectrum of B and B* is empty. 

Lemma IV.5 o'p{B) = ap{B*) = 0. 


Proof. Take any /r £ C. Then, /i £ (Jp{B) U (jp{B*) if and only if at least one of the differential equations 



:iaw sind(a;) 
V Aa; 



T[u]{x) = 0 , 


d E.k 

dx y/Ax sind(a;) 


^aw sind(a;) 



= 0 
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has a square integrable solution. Let E = 1 + a and recall that k = k + 1/2. The solutions of these differential 
equations are 


(x) = (tan 


'd{x) 




'>Pm (x) = ce 




tan 


i9(x) 


k+i 


with 

arctan (^/a cos d (a;)) 

7^ ■ 

The functions (p[^] and are defined up to a multiplicative constant c G C. Without loss of generality we set c = 1. 
Let us show that this functions are not square integrable on the interval I = ( 0 ,X 2 )- We start by observing that 

(¥>[/.], ‘PM (x)‘Pm (x) = 

with 


‘PM 


dd 

V 


aPW 


tan — 

2 


-(2fc+l) 


d-d 


f{x) = —pLX + ( ^ + :^ ) ■\/a-arctan(-\/CTCOsd{a;)) + aw(l + a) 


p{'d) = — 2 {^p)-d + { 2 k + 1)i/ct arctan {y/acos-d) + 2 au){l + cr) 


arctan {^/acos-d) 

7 ^ 


Taking into account that 1 < a/A^ < y/l + a for -d S [0,7r] and that the function arctan cos d) is continuous and 
decreasing on [0, tt] let M = inf^g[o_ 7 r]{e^^’’^} > 0. Then, 


{‘PM^^P]) > 


M 


^/T 


r '^/2 


■d 

tan- 


-(2fc+l) 


dd + 


'■ 7^/2 


d 

tan- 


-(2fe+l) 


dd 


For fc > 0 we have tan ('d/2) < 1 when -d G (0,7r/2) and therefore 


/o 


■^/2 / ^ 

tan — 

2 


-(2fe+l) 


dd > 


>0 


^ 

tan — 

2 


dd = 2 In sin ■ 


7/2 


= + 00 . 


When A: < —1 we have tan ('d/2) > 1 on (7r/2, tt) and hence 


/ir/2 


tan ■ 


d 


-(2fe+l) 


dd > 


/2 


d d 

tan — dd = —2 In cos — 
2 2 


= + 00 . 


-K {2 


In any case it results ^ L'^{I,dx). Similarly, it can be shown that '0[^] ^ L'^{I,dx). □ 


For p, S C we introduce the formal differential expression defined by 05^ = 


0 «+ - + 
«8_ - 71 0 


and we associate to *8, 


the differential operator Bfj,S = 25^5 with = D{A). Moreover, in the notation of the previous remark we have 

B — Bo = . The operator is self-adjoint for any p gC because A is self-adjoint and Bfj, — Ais symmetric 

and bounded. Following an analogous proof to that of Theorem 2.14 in it can be shown that aess{B^) = 0. This 
fact together with the next result implies that B^ is boundedly invertible. 

Lemma IV. 6 For all p G C we have 0 ^ (Tp(i3^). 


Proof. We prove it by contradiction. Suppose 0 G ap{B^) and S be an eigenfunction of Sp with eigenvalue 0. Then, 


0 = B^S = Bp 





S 


implies that either B — p or B* — p, is not injective which is in contradiction to crp{B) U ap{B*) = 0 following from 
the previous lemma. □ 
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We now derive an auxiliary result that will allow us to prove that the operator B has compact resolvent from 
which it will follow that the angular operator has compact resolvent as well. According to Lemma II V.6I we have 
0 £ C\{ap{Bf_i) U (Jess{Bp)) = p{Bf_i), where p{Bp) denotes the the resolvent set of B^. Hence, B — p, and B* —Ji are 
boundedly invertible. Moreover, their resolvents and the resolvent of Bp are connected as follows 

(B-p)-^ 0 \_fB-p 0 

0 {B*-p)-^) 0 B*-pj 


In particular, we have that the ranges oi B — p and B* — p are such that rg(i3 — p) = Tg{B* — p) = Lp‘{I, dx). Hence, 
we have shown that ij{B) = <j{B*) = 0. 

Lemma IV.7 Let p G C and y)[^] and he defined as in the proof of Lemma \IV.5\ Then, the operators (B — p)~^ 
and {B* — 'p)~^ map functions g,h G dx) to 


1 f Y 9 (t) dt ifk >0 

\ ipfp^(t)g(t) dt ifk< -1 


(43) 


and 


{B* — p) ^h{x) 


1 

V'M (a^) 


*/ fc >0 

Sx 2 '^[Y)Ki) dt ifk<-l 


X G I. 


(44) 


Proof. From the proof Lemma IIV.5I it follows that is a solution of {B — p)u = 0 and that (p[p] is a solution of 
{B* — p)v = 0. To prove (l43t and (l44)) we first show that {B — p)G{x) = g{x) and {B* — p)H{x) = h{x) with x G I 
hold formally, where G and H denote the r.h.s. of (H51) and (H^ . respectively. Let us start by assuming A: > 0. Then, 
for g £ igiB — p) we find 


{B-p)G{x) = (B-p) 




J TliJ.]i*)9{t) dt^ ={B- p) (^ij[p]{x) J p[p]{t)g{t) d?j , 


= {B'f[p]{x) - pip[p]{x)) I (pip]{t)g{t) dt + ilJip]ix)— I (p[p]{t)g{t) dt = ■f[p]{x)(pipi{x)g{x) = g{x) 


where we used (a;)(/?[^] (cc) = 1 and (f8+ — p)'fip] = 0. The case fc < —1 and the equation for h can be 

shown in a similar way. Let us prove that G G D{B) and H £ D{B*). We will do it only for G G D{B) 

in the case k > 0 because the statement for fc < —1 and H G D{B*) can be proved analogously. Since 

D{B) = {g G L'^{I,dx) I g absolutely continuous, f8+g G L‘^{I,dx)}, we need to verify that G G L‘^{I,dx). Let 

k > 0. Then, 


llGiii. = 


|G(x)p dx = f 

Jo 


\T[u] i^W 


nu](t)9i0 dt 




|5(t)I 


dr 


dd 


\J A,- j ^/A^ 


< 


dx < 




Tlu] W 


nu] 0 


\git)\dt] dx, 




| 5 (r)|dr dd 


(45) 


with 0 < d < T < TT. Note that by assumption g G L^((0, tt), dd/ A^) and its restriction on the interval (0,i?) will 
belong to L^((0, d), dr/v^A^) for d G (0, tt). Furthermore, we have 




2 



2fc+l 


with 


x(t, d) = 2 



[arctan(-\/o'COsr) — arctan(-\/CTCOsd)] . 


First of all observe that e with K = max{l, Furthermore, the 

function x(t, d) has an absolute maximum on the square [0, tt] x [0, tt] at the point (0, tt) where x(0, tt) = 2 arctan-ya. 
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Hence, there exists a positive constant C such that < C with C = ^)arctanVCT^ Finally, since 

0<i?<T<7rwe can show that 


tan f d 

- - < - 

tan ^ T 

as follows. Since the above inequality is equivalent to the inequality ^ ^ , we consider the function / : 

tan. — 

(0,7r) —>• R such that f{'d) = ^ ^ . Clearly, / is continuously differentiable and the inequality we need to prove is 

equivalent to /(i?) < /(r) for 0 < d < r < tt. Hence, it suffices to show that / is a monotonously increasing function 
of d. A straightforward computation shows that 


/ w = 


1 ( sinf A 

2dcos2|l I 2 


since (sin|)/(d/2) < 1 and cos ^ < 1. Finally, d < r implies that 


tan ^ 
^ < 1 
tan 77 


and we conclude that 


{t) 




(46) 


Hence, yj[^](r)/y3[^](d) £ L‘^{{0,TT),dT) for each fixed d G (0,7r). Therefore, we can use the Cauchy-Schwarz inequality 
to estimate the inner integral in (l45l) as follows 


f r 


|g(T)|dr < 



ntA (t) 


< CK f f 

J-d J'd Jo 


|5(r)|^dr < C'o||g||i2e^l*^'"l’^(7r-d) 


with Co = CK. Inserting the above expression into (|^ yields ||G |||2 < ACo/ 2 )e^'^^^'^'^\\g\W '2 < +oo and this 
concludes the proof. □ 


Note that Lemma FlV.?! gives explicit expressions for the inverses of B and B* by choosing /i = 0. Now that we have 
obtained an explicit form of we can show that B, and therefore the angular operator A, has compact resolvent. 

Lemma IV.8 The operator B has compact resolvent. 


Proof. Tho show that the operator B ^ 


0 

H-i 0 


is compact, it suffices to show that the operators B ^ 


and {B*)~^ are compact. We limit us to prove that B~^ is compact in the case k > 0 because the case fc < —1 and 
the corresponding assertion regarding {B*)~^ can be proved in an analogous way. From the previous lemma we have 
for fc > 0 and g £ L^(/, dx) 


1 

B~^9{x) = - ^ (p[o]{t) 9 {t)dt, xGl. 

T[o][x) Jx 

For each n G N we define the operators 


r„ : L\I, dx) L\I, dx), Tr^fix) = 


r+i ^[0] W/(i) dt if X G [^,X2 - i] 


V[0] 


and 


r„ : 


1 1 ' 

-,X2 - 

n n 


dx —> K 


1 1 

-,X2 - 

n n 


■‘‘"I’ = 
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These operators are bounded for all n £ N and the operators Tn are compact since the integral kernel is continuous 
and bounded (see example 4.1, Ch. Ill in [l^). For every / G L‘^{I,dx) the restriction / of / on the compact 
interval [^,X2 — lies in ([^,0:2 — It is clear that for any convergent sequence {gn)n£N with terms in 

([i,X2 — ,dx) also the sequence (g^rt)ngN with terms in ([^>^2 — ,dx) converges where 

~ . Jo if 2:^ 1 ^X 2 -^] 

9 n{x) = i , , „ n 

I 9{x) ifxG [-,X2--] 

for all n G N. Let (/m)7neN be a bounded sequence in L‘^{I,dx). Then, (/m)meN is also a bounded sequence in 
([^,^2 — ,dx^. Hence, for any n G N the sequence {Tnfm.)m&n contains a convergent subsequence. Therefore, 

also (Tnfm)men contains a convergent subsequence since Tnfm = Tnfm- This shows that the operators T„ are also 
compact. The proof is completed once we have shown that T„ —>■ B~^ as n 00 in the operator norm, that is we 
have to verify that ||T„ — B~^\\ —^ 0 as n —>■ 00. To see that, we note that for all / G dx) we have 

||(T„ - B-^)f\\l. = r |(r„ - B-^)f{x)fdx = /" \B-^f{x)fdx+ 

Jo Jo 


\T^f{x)-B-^f{x)\^dx + 


rX2 |•X2 

2 — -/ I D —1 J’/'.„M2 I / I D —1 


\B ^f{x)\^dx= / \B ^f{x)\^dx + 
Jo 


B f{x)fdx 


using the definition of Tn and B Taking into account that for all (a, h) C (0, X 2 ) we have 


f\B-^f{x)\^dx= f 

J a J a 


<P[0](a:) 


no] 


dx < 


'a \J X 


no] ft) 

no] 


\f{t)\dt dx 


</::(/; 

where we used the fact that (/3[o] is real-valued. Employing the Cauchy-Schwarz inequality and (I46|) yields 

rm 




‘^dT]d-d<CK I [I IdT 

i9(o) \J^ 


|/(r)| dr dd 


<CK I dr / |/(r)pdr dd = rCiF(d(6)-d(a))||/||i. 

Ji}{a) \Jo J \Jo J 

Hence, 

where we used the fact that d(0) = 0. Finally, 


d -I- d(x2) - d 

and taking into account d = d(x) is continuous and therefore sequentially continuous, we find that 

lim ||T„ - B~^\\ < \/nCK [d(0) -I- d{x 2 ) - ■d(a; 2 )] = 0 



Tn-B-H =< JnCK 


d ( — j -I- d{x2) — d (x2 - 


\h 


This completes the proof. □ 

Theorem IV.9 The angular operator A has eompact resolvent. 














20 


Proof. We know that both A and B are self-adjoint and therefore their spectra are real. Let us take any fi G p(A) 
and ly € p{B), then the second resolvent equation 

{A — p)~^ — {B — v)~^ = {A — p)~^{B — A + p — v){B — v)~^ 

but {B — v)~^ is compact and {A — p)~^ and {B — A + p — v) are bounded, hence the operator on the r.h.s. of the 
above expression is compact. This implies that also {A — p)~^ must be compact. □ 

We conclude this section by deriving a set of necessary and sufficient conditions for the angular eigenvalue problem 
to have polynomial solutions on the interval (0, tt). The same conditions can also be used to compute the corresponding 
eigenvalues. To this purpose we start with (l37l) and make the substitution 5'_(d) = (sin'd)“^/^S'_(d) leading to an 
ODE of the same form as dazi) with S- and C± replaced by S- and 


C± = — ± — auj sin i? 

M I 


sini? 


respectively. In the case a = A = 0 the equation can be solved and the solutions can be expressed in terms of Jacobi 
polynomials [3l|, [s^ . If A = 0 the same equation can be transformed into a generalized Heun equation [3l| . Note 
that the initial choice of working with Carter’s tetrad leads to expressions for the operators C± that are much simpler 
than the corresponding ones for the operators £ 1/2 and appearing in [s^ where Kinnersley’s tetrad was adopted. 
Finally, this ODE becomes 

S- = 0, 

where a prime denotes differentiation with respect to d and 

Hid) = —^ ( aw sin-d ) • 

^/Ks y Sind j 


A ^ 


auie^i) sind 
A + arrif. cos d 


2A,, 




A ,5 


arUe^H sind 
A + arrie cos d 


-I- A^ — a^rue cos^ d 


Introducing the variable transformation a; = (1 J- cosd)/2 which maps the interval ( 0 , 7 r) to the interval (0,1) the 
above differential equation becomes 


<PS. 

dx^ 


X — 1/2 Ax 
x(l — x)~^ 2Ax 


2 ame 

dS- 


^ 2kx — k 

2ameX J- A — arue 

dx 

4x2(1 - x)2 

A,, 


Sfc" 


auj^Ci — 2x) /S.X ^oru) 

— auj—^ — 


2 , , 2—2 


2 ame 


2 aL 0 


c)! - x)Ax 


Al Al 


2ameX -I- A — anie 


2x(l — x)Ax 


A^ — a^TOg(2x — 1)2 

x(l - x)Ax 


S- = 0, 


where a dot denotes differentiation with respect to x and A^, = 4(S — l)x{x — 1) J- Let Q denote the variable 
coefficient multiplying S- in the above equation. In order to kill terms like 1/x^ and 1/(1 — x)^ appearing in the 
expression for Q we observe that 


4x2(1 — x)2 y A 


( 2kx — k 


Al ^ Bi B2 ap + aix 

X x2 1 — X (1 — x)2 Ax A2 


with 


Al =- 

2S 


fc(5S-4) -h2(S- 1) 


, k{k + l) ^ k 
A2 = - - -, B, = - — 


A:(5S-4)-2(S- 1) 


= - 


k{k-l) 


ao = -^(2fc + l)(S-l)^ = C, = -4p{E-ir 


and we rewrite Q as follows 


, A2 B2 Ax[Ai + {Bi - Ai)x\ + x{l - x)[ap + aix) A aioE{l - 2 x) 


x(l - x)Ax 
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C\ — 

a! 


2 ame 


2 auj 


2ameX + A — amg 


2 x{l — x)A^ 


+ 


— a?m 1 { 2 x — 1)^ 
x(l - x)A^ 


At this point we introduce the s-homotopic transformation S-{x) = x“(l — x)^^{x) with a,/3 G C and we find that 
the equation satisfied by $ is 

TT + P{x)^ + q{x)<^ = 0 

ax^ ax 

with 


p{x) = 


2 a + l/ 2 - (l + 2 a + 2/3)a 
a;(l — x) 


2 amg 


A, 


2Ax 2ameX + A — amg ’ 


q{x) = 


-a/2 P^-/3/2 2aP + a/2 + /3/2 aA^ 


/3A, 


+ Qix). 


^ 2ame[(a +/3)x — a] 

x"^ {1 — x)"^ x{\ — x) 2xAx 2(1 — a;)Aa; a;(l — x)(2ameX + A — arue) 

The coefficients of x~‘^ and (1 — x)~‘^ will vanish whenever 

,_a_t(t + l) ^ g_i:(t-l)^ 

2 4 2 4 

The only acceptable roots are those keeping the solution regular at x = 0 and x = 1 and they are given as follows 

1 1 


“"4 + 2 




l + /3=i + i 

4 ^ 2 ’ ^ 4^2 


k - 
2 


i , M 

4 2 ' 


Finally, multiplying the equation for $ by x(l — x)A^(2ameX + A — anig) yields 

(72 (j) /7(j) 

X{x)^+Y{x)^ + Z{x)<^ = 0, 
dx^ dx 

where X, Y, and Z are polynomials of degree 7, 6 , and 5, respectively, given by 


(47) 


A(x) = x(l — x)A‘^(2amgX + A — arug) = ajX^ 

j=o 


with 


07 = —32ame(S — 1)^, oe = —16(S — 1)^(A — lanig), 05 = —16ameS(S — 1) + 48(S — 1)^(A — Sanig), 

04 = —8S(S — 1)(A — Some) — 16(S— 1)^(3A — 5amg), 03 = —2amgE? + 16S(S — 1)(A— 2amg) + 16(S — 1)^(A — amg), 
02 = —S^(A — Same) — 8S(S — 1)(A — ame), oi = S^(A — ame), oq = 0, 


6 


Y{x) = X{x)p{x) = ^ bjX^ 
j=o 


with 


&6 = -S2amg{E - 1)2(1 + 2a + 2/3), 65 = 32(S - l)^[amg{7a + 5/3 + 3) - A(1 + a + /3)], 

64 = — 8 (S — 1) {(S — l)[5ame(8a + 4/3 + 3) — 2A(6a + 4/3 + 5)] + ame{4a + 4/3 + 1)} , 


63 = 4(S - 1) {(S - l)[20ame(3a + /3 + 1) - A(28a + 12/3 + 19)] + 4ame(5a + 3/3 + 1) - A(4a + 4/3 + 3)} 
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62 = —4ame(a + /3) + (^ — 1) {S[—10ame(10a + 2/3 + 3) + 2A(32a + 8/3 + 17)] + 4ame(7a — /3 + 4) — 16A(1 + 2a)} , 

&i = 2aTOe(3a+/3)—A(2a+2/3+l)+(S— 1 ) {S[2ame(lla + /3 + 3) — A(18a + 2/3 + 7)] + 2ame(3a + /3) — A(2a + 2/3 + 1)} , 

"2 

60 = + 1 )(A - aWe), 

and 

5 

Z(x) = X(x)q(x) = CjX-^ 

3=0 


with 


C 5 = — 16aTOe(S — 1)[(S — l)(k^ + a + /3 + 4a/3) + 2a?m1], 

C 4 = 8 (S — 1 ) |(S — l)[fc^(5ame — A) + ame{2k + 2aoj + 5a + 5/3 + 20a/3) 
—A(4a^ + 3a + 3/3)] + 2a^me{7>am\ — Xrrie + w)} , 


C 3 = Sa^rrieiuJ^ — ml) + (^ — 1) |^[8A(2/c^ + k + 8a/3 + 7a + 5/3) — 8 ame{ 6 k^ + 4fc + 20a/3 + 5a + 5/3)+ 

8 a'^me 0 j{au! — 4)] — 8\{2k‘^ + k + 8a/3 + 7a + 5/3) + 8ame(5fc^ + 4fc + 16a/3 + 5a + 5/3 + A^)+ 

8a^me(4A77ie + — lOamg)} , 


C 2 = 4a^(A — 3ame)(co^ — m^) + (S — 1) |s[—4A(4fc^ + 3fc + 12a/3 + 12a + 6/3) + AameiSk^ + 5k + 20a/3 + 5a + 5/3)+ 
4a^a;(4me + Aw — 3amew)] + 4A(3A:^ + 3fc + A^ + 8a/3 + 10a + 4/3) — 4amg{5k‘^ + 5fc + 8a^ + 6a + 4/3 + 3A^)+ 

4a^[A(a;^ — 6mg) — ame{3uj^ — lOWg)]} , 

Cl = —4a^A(w^ — TOg) — 2a^mg(3mg — 2w^) — ameik"^ + 4a/3 — a — /3 — 2A^) — 2Aaa;+ 

(S — 1) |s[2A(4fc^ + 3fc + 8a/3 + 9a + 3/3) — ame(13fc^ + 4fc + 20a/3 + 5a + 5/3) + 2aw(2a^a;mg — 2awA — A)] 

—4A(fc^ + fc + A^ + 2a) + aTOe(3fc^ + 4fc — 4a/3 + 9a + /3 + 6A^) — 2a(wA + 5a^m\ — 2a^a;^TOg — 4aXm1 + 2aAa;^)| , 


Co = — |(S — 1)[—A(5A:^ + 2k + 4a^ + 5a + ^) + ame{5k^ + 4a/3 + a + /3) + 2aa;(A — am-g)] 

—X{k‘^ — 2X^ + 4a/3 + a + /3) + ame{k^ — 2k + 4a^ — 3a + /3 — 2A^) + 2a(A — ame){uj — aTOg)| . 

Note that the coefficients 07 , • • • ,ao are not algebraically independent since their sum vanishes. 

Theorem IV. 10 If a = 0, equation J^Tj) admits a constant polynomial solution. 
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Proof. Without loss of generality let $( 0 :) = 1. This will be a solution of (|Tfl) whenever Z{x) = 0, that is C 5 = • • • = 
Cq = 0. From the expression giving C 5 we see that it will vanish if a = 0 for which S = 1 or rrie = 0. In the case a = 0 
all coefficients Ci vanish identically and therefore $(x) = 1 is a solution of (iTfl) . If me = 0 and k > 1/2, we find that 
C 5 = —4(S — I)^(2fc + 3)(2A; + 1) can never vanish. The same situation occurs for k < —1/2. □ 

Remark IV. 11 The case 0 = 0 corresponds to the Schwarzschild-deSitter metric for which one can show that 0 
becomes a hypergeometric equation and the eigenvalues can be computed explicitely as given in ^ 1 , [H/. 

To find a set of necessary and sufficient conditions for the existence of polynomial solutions of the form $(x) = x —xi, 
we generalize the so-called functional (or analytic) Bethe Ansatz method used in so that it can be applied to 
equation 

Theorem IV.12 The polynomial <I>(x) = x — xi is a solution of if and only if 

be +C 5 = 0, 

bexi -I- &5 -I- C 4 = 0 , 

bexl + & 5 X 1 -I- 64 + C 3 = 0, 

bexl + b^x\ + 64X1 -b 63 + C2 = 0, 

be,x\ -b h^x\ + bix\ -b b^xi -b &2 + ci = 0 , 

where xi is determined by the equation Y(xi) = 0. 

Proof. First of all, by replacing the ansatz <i>(x) = x — xi into (iTTll we find that Y{x) + Z{x){x — xi) =0 and hence 

6 5 

—Co = (x — Xi)“^ biX^ -b CiXh 
2—0 2=0 

The l.h.s. is a constant while the r.h.s. is a meromorphic function with simple pole at x = x^ and a singularity at 
X = 00 . The residue of —cq at x = Xi is 

6 

Res(—Co,X = xi) = lim (x — xi)(—Co) = 7 bix\=Y{xi). 

X—^X\ • ^ 

2=0 


Then, we have 


6 5 

—Co — {x — xi)“^Res(—Co, x = xi) = {x — hi{x‘^ — x\) + c^x* 

2=0 2=0 

and after simplification of x — xi in the denominator of the term appearing on the r.h.s. of the above expression we 
end up with 


with 


5 

- Co = ^ liX^ + 
2=1 


Yjxi) 

X — Xi 




4 3 2 

='^bi+2X\+Ci, l2='^bi+3,x\+ C2, ='^b^+ix\ + c^, 

2=0 2 = 0 2=0 


(48) 


4 

74 = h^Xi + ^5 + C 4 , 75 = ^6 + C 5 , K = h^ix\. 

2=0 

The r.h.s. of (H51) is a constant if and only if the 7 /s as well as F(xi) vanish. Note that when this set of equations 
are satisfied, then —co=K. □ 
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The same strategy can be applied to show that <i)(a;) = (x — Xi){x — X2) will be a solution of (iTTl) whenever the 
following set of conditions are satisfied, namely 

2(07 + feg) T C5 = 0, (2a7 + bfj'j^Xi + X2) + 2(05 + 65) + C4 = 0, 


2(07 + a 5 ){xl + X1X2 + xl) + (2a6 + b 5 ){xi + X2) + be{xl + xl) + 2(05 + 64) + C3 = 0, 

(ll + 3^2) [207(^3 “f ^2) ^4] “f 206(3^4 + X\X 2 + x^) + bfy^Xi + X2) + b^(xi + X^') + 2(04 + 63) + C2 = 0 , 

2 ar{xf + xfx2 + xfxl + Xixl + X2) + ( 3 ;i + 3 ; 2 ) [206(3:4 + X2) + 204 + 63] + bQ{x\ + 3:2)+ 

65(3:4 + 3:2) + 64(3:4 + 3:2) + 2(03 + 62) + C4 = 0, 

207(3:4 + x\x2 + x\x\ + 3:43:2 + 3:43:2 + 3:2) + 2 ae(xf + xfx2 + xfxl + XiX^ + 3:2) + (3:4 + 3:2)[2a5(3:4 + xl) + 203 + 62] + 

204(3:4 + 3:43:2 + 3:2) + 65(3:4 + 3:2) + 65(3:4 + 3:2) + 64(3:? + 3:2) + 63(3:4 + 3:2) + 2(02 + 64) + Co = 0 

together with Res(—co,3: = Xn) = 0 for n = 1, 2. This equation gives rise to two algebraic equations for the roots Xn, 
more precisely 

^ _l_ _ Q 

m = l ^(3:71) 

m^n 


V. THE RADIAL SYSTEM 


In this section we show that the Dirac equation in the non-extreme Kerr-deSitter metric does not allow for bound 
state solutions by computing the deficiency index of the radial operator associated to the system (1551) . More precisely, 
since the deficiency index counts the number of square integrable solutions, it suffices to show that the deficiency index 
of the radial operator vanishes. We start by observing that the components of the radial spinor satisfy the relations 
R- = i?+ and i?+ = R- as it can be easily verified from (j35p . This property motivates the ansatz (r) = F{r)—iG{r) 
and R+{r) = F{r) + iG{r) used in the proof of the following theorem. 

Theorem V.l In the non-extreme Kerr-deSitter metric the Dirac equation does not possess bound state solutions. 

Proof. The proof relies on an application of Theorem 5.2 in and follows the strategy adopted in the proof of 
Theorem 5.1 in [4l|. To this purpose we cast the radial system (l3^ into the form 




■ V(r 


dr ^ " VA7 

imeX — A 


—iruer — A 


V Ar 


d 

dr 



aYIlI 


= 0, V{r)=E uj{r‘^ + + ak 


(49) 


In order to bring (1491) into a Dirac system we transform the dependent variable according to R-(r) = F{r) — iG{r) 
and R+{r) = F(r) + iG{r) and we introduce a new independent variable u defined through the relation 


Then, (l4^ becomes 


du E(r^ + a^) 
dr Ar 


(iim(u) := -f B(u)n = ujn 
du 


(50) 


(51) 



" 0 lA 
-1 0 j 


/ meTy/ A.r — a'E.k 

\ 

^=( 

. B{u) = j 

AVA7 

H(r^+a^) _ I 

meT\/Ar+aSfe / 







with D = {F, G)^ and 
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where the radial variable r is now a function of the tortoise coordinate u. In the non-extreme case the behaviour of 
the solution of (1501) in proximity of the event and cosmological horizons is captured by the following formulae 


'^(^) = -In. (r — r+) -I- 0{r — r+), 

2av_|_ 


nC?”) = X— In [tc -r) + 0(r - Tc), 

ZKn 


(52) 


where and Kc denote the surface gravity at the event and cosmological horizon, respectively, and are given by 

3(r+ — M) — Ar+(2r^ -|- a^) Arc{2r‘^ + a^) — 3(rc — M) 

^ 3S(r2 + a2) ’ """ ^ 3S(r2 + a^) ' 

Since > 0 on the interval (r_|_,rc), it follows from (1501) that u must be an increasing function of r and therefore, 
the numerator in the expression of is positive. This means that u —>■ —oo as r —>■ r^. By the same token we 
can also conclude that the numerator in the expression of Kc must be negative and therefore, u —>■ -l-oo as r —>■ rc- 
Furthermore, u maps the interval (r_|_,rc) to the real line. Taking into account that the radial spinor R = 
is square integrable on (r+,rc) if 


{R,R) 



■\R? 



-I- a? 
Ar 


R-R+ < -1-00, 


it results that the transformed spinor is square integrable whenever 


{n,n) 



du [F{u)‘^ + G'(u)^] < -boo. 


(53) 


The formal differential operator il is formally symmetric since J = —J* and B = B* where star denotes complex 
transposition. Let Smin be the minimal operator associated to il such that Smin acts on the Hilbert space du)^ 

equipped with the inner-product (l53)) . Then, Smin with domain of definition D{Smin) = (^^(K)^ such that Smin^t := 
ilH for H G D{Smin) is densely defined and closable. We denote by S the closure of Smin and apply Neumark’s 
decomposition method [d^. To this purpose, let Smin,± be the minimal operators associated to il when restricted 
on the half-lines [0,-|-oo) and (—oo,0], respectively. We consider Smin,± on L‘^{R±,du)'^ equipped with (l53l) . The 
operators Smin,± with domain of definition D{Smin,±) = (^“(K-i)^ and Smin,±^± ■= ilfl± for G D{Smin,±) are 
densely defined and closable. Furthermore, il is in the limit point case at ±oo. This can be seen as follows. First of 
all, we recall that since the limit point and limit circle cases are mutually exclusive, we can determine the appropriate 
case if we examine the solution of (ISTj) for a single value of w. Hence, withont loss of generality we set w = 0 and 
consider the system 

dfl , ^ ^ 

J— + B{u)n = 0. (54) 

du 

First of all, we observe that the matrix B[u) converges for u —>■ —oo to the constant matrix 


Bq = lim B{u) 

lA—>• —OO 


ak 

r^+a^ 

0 


0 

ak 

7^ 


This observation together with (1521) suggests an asymptotic expansion of B{u) in powers of 6"^+“. A straightforward 
computation gives 


B{u)=Bo + B{u), \B{u)\<Ce^+'^ 

for some positive constant C. Finally, applying Theorem 1 in Ch. IV in [d^ yields that the system ((5^ has asymptotic 
solutions for u ^ —oo given by 


F{u) = cos 


ak 


rl + 


{l)+0 (e”*“). Glu) = sin ( ; ) + O (e-“). 


The problem at the cosmological horizon can be treated similarly and for u —>■ -l-oo we find that 


F{u) = cos 


ak 




J ^ +0(e'=‘=“). 


(55) 


(56) 
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By inspecting (1551) and (1561) we can immediately conclude that the differential operator it is in the limit point case 
at ±c». This implies that the operators Sm,in,± are essentially self-adjoint. Let S± denote the closure of Smin,± and 
N±{S±) the corresponding deficiency indices. If v± denotes the number of positive and negative eigenvalues of the 
matrix iJ, then, since = 1 = v-, Theorem 5.2 in implies that N±{S+) = 1 = N±{S-). Since zero is the only 
solution of (iMl) in du) the orginal system (ISTI) is definite on R+ and ]R_ in the sense of [i^. Finally, (5.11a) in 
Proposition 5.4 in yields that the deficiency indices for S are 

N±{S) = N±{S+) + N±{S.) -2 = 0. 

This implies that the radial system m does not possess any square integrable solution on the whole real line and 
this completes the proof. □ 


Appendix A: Derivation of formulae (1131) and (1141) 

First of all, the quartic equation ([2]) is already in reduced form. Let with i = 1, • • • ,4 denote the roots of ([2]). 
Applying Vieta’s formulae we find the following useful relations among the roots of A^, namely 

ri+r2+r3+r4 = 0, ri(r2 + + r4) + r2(r3 -|- r4) -I- r3r4 = p, 


fir2r3r4 = u, ri (r2r3 -I- r2r4 -I- r3r4) -|- r2r3r4 = -q 
with p, q, u defined in ([21). Let us consider the following polynomials in r^, namely 

^1 = (ri-I-r 2 )(r 3 -I-r4), ^2 = (ri-|-r 3 )(r 2 -b r4), = {ri + r4){r2 + r^). 

Using Vieta’s formulae we find that 

zi + Z2 + Z3 = 2(rir2 -b rir3 -b rir4 -b r2r3 -b r2r4 -b r3r4) = 2p, (Al) 

Z\Z2 + Z2Z3 ^ Z3Z1 = p"^ - Au, ziZ2Z3 = -q^. (A2) 

This means that zi, Z 2 , and Z 3 are roots of the so-called resolvent cubic z^ — 2pz^ + {p^ — Au)z -b = 0. Since ([U is 
reduced, Vieta’s formula ri+r 2 + r 3 + r 4 = 0 holds and equations and (IA2I) can be solved yielding ri -br2 = \/—z\, 
r 3 + r 4 = -yj-zi, n -b r3 = yJ-Z 2 , r 2 +r 4 = -^-Z 2 , ri -b r4 = y/-Z 3 , and r 2 + rs = -^-zj . The ambiguity in 

the choice of the sign of the square root can be fixed according to the relation yj—z\^—Z 2 \/—Z 3 = —q. Taking into 

account that in the present case 9 > 0 we choose the sign of the square root so that 

2ri = -yJ-Zx - yJ-Z2 - \/-Z3, 2r2 = -yJ-Zx -b yJ-Z2 -b \/-Z3, 

2r3 = ^yj-zx — \/-Z2 + 'J-Z3, 2r4 = +\/-Zx + \/-Z2 - yJ-Z3. 

To find the roots zi, Z2, and Z3 of the resolvent cubic we reduce it to the special cubic -b pw -b g = 0 by means of 
the Tschirnhaus transformation z = w -\- (2p)/3 where 

/P" . ^ ^ 2 3 8 2 

P=-ly+4Ml, + 9 • 

Let wx, W2, and W3 denote the roots of the special cubic. With the help of Vieta’s formulae 

Wx + W2 + W3 = 0, WxW2 + W2W3 + W3WX = P, WXW2W3 = —q 
the discriminant of the special cubic turns to be 

D = {wj — wiY = — 4 :{wxW 2 + VJ 2 W 3 -b W 3 Wx)^ — 27{wxW 2 W 3 )^ = —4^ — 27q^. 

l<j<i<3 

Let uj = be the primitive root of unity and introduce the Lagrange substitutions 

1 ^ 

a= fc = 0,l,2. 
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Using = 1 and wi + W 2 + W 3 = 0 we find that 

3^0 = Wl + W 2 + W 3 = 0, 3^1 = Wi + UJW2 + Uj'^Ws, 3^2 = Wi + 0j'^W2 + UJW3. 

A somewhat tedious computation gives 



where we used the fact that + w + 1 = 0. Hence, we obtain 



where the signs of the cubic roots must be chosen so that ^ 1^2 = —p/S. The Lagrange substitutions introduced earlier 
can be solved uniquely in terms of and ^2 and we obtain 

2 

zc, j = 1,2,3. 

fc=i 


Finally, we get the roots of the special cubic as 

Wl=Cl+?2, W2 = + U}^2, W3 = LU^l + U}^^2- 


Appendix B: The extreme case 

Let ^ and factorize the polynomial in (jd]) as (p — p _ ){p — p-){p — pcY which compared with (jd]) gives rise 

to the following system of equations for the unknowns p _, p_, and pc 

p — + p- + 2pc = 0, p — p- + 2{p — + p-)pc + Pc = — 1, 


-2p__p_pc - (p— + P-)Pc = Pi.+, P-P-pI = -c?- 

Note that since p _< 0 and > 0, the last equation implies that p_ > 0. Moreover, the first equation permits to 

express p _in terms of p_ and pc so that the above system reduces to 


P- (p- + ‘^Pc) + 3Pc 

= 1 — 

(Bl) 

2p-Pc(P- + ‘2‘Pc) + ‘^pI 

= Pi,+ , 

(B2) 

P-Pc(P- +2pc) 

= a". 

(B3) 


We have to deal with an overdetermined system and special care must be exercised in analyzing its solutions. First 
of all, we note that all three equations above are quadratic polynomials in p_. Solving (IbT])-(Ib 31) with respect to p_ 
and taking into account that p- > 0 , we find 


5-,/ = -pc + -2p\ 


P-,ii = 


~2Pc + \/2pi,+Pc 

2.pc 


P-jii = 


\fpt 




At this point we must require that + 2p^ < 1. This condition is always satisfied because > 0 and therefore 
\/l — — 2p1 > Pc- Squaring we get 1 — > 3pl and hence 1 — — 2pl > 2,pl — 2p1 = > 0. Furthermore, the 

Roman numeral attached to the roots indicates that p-j is a root of m and so on. In order to make the system 
(iBB-dBl consistent, we must require that p_ j = p-ji = p-jij. This implies that 



(Bd) 


x/l - - 2p2 = 
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From the first equality in (IB4I) we obtain the biquadratic equation 

3pt - (1 - a^)pl + a^ = 0. (B5) 

Since pc > 0, the only two acceptable solutions are 


Pc,± = - 6a2 ± (a2 _ a'i){a‘^ - a%) 


with a± = 2 ± Equating the first term and last term in (jB4l) and then equating the second and last term in (IB4I) 
we find that 

Mi,+ = 2 : m.+ = 2(1 - a^)pc - ^pl- 

We see immediately that the above expressions for pi + are equivalent provided that pc is a root of (IB5I1 . Finally, 
a lengthy but straightforward computation shows that pi.+ coincides with (??) only if we choose pc = Pc,+ - The 
solution Pc = Pc,- must be disregarded because would lead to the contradiction /ii_+ = p,2,+ - Concerning the case 
P- = p+ and /i = p2,+ we observe that the corresponding set of equations involving p+ and pc can be obtained 
directly from (iBT])-(imi) with p_, Pc and replaced by pc, p+, and p2,+, respectively. The rest of the proof is too 
similar to the previous case to be presented here. 
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